New magic number for neutron rich Sn isotopes 
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The variation of E(2~t) of 134-140 5n calculated with empirical SMPN interaction has striking 
similarity with that of experimental E(2f) of even-even ls ~ 22 and 42_48 Ca, showing clearly that 
N=84-88 spectra exhibit the effect of gradual filling up of u(2ff/ 2 ) orbital which finally culminates 
in a new shell closure at N=90. Realistic two-body interaction CWG does not show this feature. 
Spin-tensor decomposition of SMPN and CWG interactions and variation of their components with 
valence neutron number reveals that the origin of the shell closure at 140 grilles in the three body 
effects. Calculations with CWG3, which is obtained by including a simple three-body monopole 
term in the CWG interaction, predict decreasing E(2f) for 134-138 Sn and a shell closure at 140 Sn. 

PACS numbers: 21.60.Cs,21.30.Fe,23.20.Lv,27.60.+j 



The evolution of shell structure away from stability [l| 
has been a topic of intense theoretical @, Q and experi- 
mental [|]-[6| studies since last decade. Theoretical stud- 
ies have identified different reasons for the phenomenon 
of shell evolution with neutron or proton excess. Among 
the various components of the nucleon - nucleon interac- 
tion, the spin-orbit, tensor part and three - body effect 
play important roles in the shell evolutions Due to 
tensor interactions, nuclear mean field undergoes varia- 
tions with neutron excess. This leads to monopole mi- 
gration Q . It is observed for both proton-rich as well as 
neutron-rich nuclei. While approaching the neutron drip 
line, the neutron density becomes very diffused Q which 
can also lead to shell quenching. For exotic light nuclei 
the well established magic numbers for the stable nu- 
clei are found to be modified or new magic numbers have 
evolved. At least four doubly magic oxygen isotopes have 



been observed 
Brown and Richtcr 



They are 14 (9, m O, 22 and 24 0. 
have framed a generalised new rule 
for magic numbers, valid specially for lighter nuclei. For 
heavier nuclei, experimental production of neutron rich 
isotopes are more difficult. There are severe limitations 
in acquiring spectroscopic information on them due to 
their low production rates and lifetimes. However, the 
shell evolution for neutron -rich nuclei above the doubly 
magic 132 Sn core is recently a topic of great interest. The 
Sn is otop es in particular, pose many interesting problems 
@> lij-lirj] in the study of evolution of nuclear structure 
with increasing neutron number. The near constancy of 
the first 2 + energy of Sn isotopes for A=102 to 130 at ~ 
1.2 MeV is a text book example for seniority conserved 
spectra [l(|. But the two valence neutron isotope of Sn 
just above 132 S' 



n, i.e.. 



134 



Sn, shows a sudden depression 
in 2f energy to 726 keV. This depressed energy is not 
only interesting from the point of view of nuclear struc- 
ture [17j . it should also have an important implication 
for the r-process scenario fl5l |. 

Large basis untruncated shell model (SM) calcula- 
tions have been done in the valence space consisting of 



7r(lff7/2,2G?5/ 2 , 2d 3/2 , 3si/ 2 , I/J11/2) and v(lh Q / 2 , 2/ 7/2 , 
2/5/21 3^3/2,3^1/2, 1*13/2) orbitals above the 132 Sn core 
usingboth realistic CWG @,|l3,[ll and empirical SMPN 
Tol . IT} (1+2) - body Hamiltonians. The Hamiltonians 
have the same set of single-particle energies [l(| of the 
valence orbitals but different sets of two-body interac- 
tion matrix elements (tbmes). It should be noted that 
the neutron-neutron part of the CWG tbmes is the same 
as that used by Kartamyshev et al. @. The shell model 
codes OXBASH and NUSHELL@MSU have been used 
[l8l |. The two theoretical results differ dramatically for 
136 ' 138 SVi pjj. The realistic interaction CWG predicts 
nearly constant energies of 2^ states for the even-even 
Sn isotopes above the doubly magic 132 Sn core, normally 
expected for singly-magic nuclei. On the other hand, the 
empirical interaction SMPN predicts decreasing E(2* ) 
energies with increasing neutron number. The calculated 
energies with SMPN fit in the systematics [l(| for the 
experimental E(2^) energies of their isotones with Z>50. 
They also agree with the trend shown by the Casten - 
Sherrill systematics for E(2f ) energy differences of Sn 
and Te isotopes having same neutron number [ToL HH [T^ . 
It has been shown in Ref. [l(| that this non-constancy of 
E(2^~ ) in Sn isotopes above 132 Sn is a strong possibil- 
ity. The prediction for dramatic decrease of the E(2f ) 
of neutron-rich Sn with increasing neutron number for 
N=84-88 using SMPN interaction was considered 10] to 
be an effect showing weakening of the Z=50 shell gap. 
But the new result, which we report in this letter, on 
140 Sn, its high 2^~ energy and its comparison with exam- 
ples from other neutron - rich domains clearly show that 
N=84-88 spectra with SMPN show the effect of gradual 
filling up of ^(2/7/2) orbital which finally culminates in 
a new shell closure at N=90. We show that the realistic 
CWG predicts similar results, that is decreasing 2+ en- 
ergies and a shell closure at 140 Sn if three body effects 
are included in it. 

The results for the E(2+) energies of isotopes of Sn 
for A=134-140 have been shown in Fig.l as a function 
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FIG. 1: Comparison of E(2~t) energies as function of valence 
neutron number for oxygen(O), calcium (Ca) and Sn isotopes. 



of valence neutrons above 132 Sn core. The experimental 
energies for 13213i Sn are shown in the figure. The pre- 
dicted energies using CWG and SMPN interactions are 
compared. With CWG interaction, as mentioned above, 
the 0^ - 2f spacing remains nearly constant at around 
750 keV for 136 ~ 14:2 Sn, except for a small increase at 
140 Sn [t| due to the filling of the v(2f 7 / 2 ) single-particle 
orbit. It has been identified Q to be a weak shell closure 
at 140 Sn. The seniority conserved character of the low 
lying spectra [l6| observed for Sn isotopes below 132 Sn 
is also preserved for isotopes above N=82 shell closure. 
The nearly constant E(2^~) value at around 1.2 MeV for 
52 < TV < 80 reduces to around 750 keV above N=82. 
The same figure also contains variation of the experimen- 
tal E(2j ) energies of even 2 oCa isotopes from A=40 to 
48 as a function of valence neutrons above 40 Ca core. 
Similarly, the variation of the experimental E(2^~) ener- 
gies of even gO isotopes from A=16 to 24 as a function of 
valence neutron numbers above 16 O core are also shown. 
The variations of experimental E(2+) with the valence 
neutron number for two different mass regions and shells, 
show striking similarity with the theoretical predictions 
with SMPN in the 132 Sn region. The gradual filling of the 
1/(2/7/2) orbital by neutrons is very distinctively shown 
by the variation of E(2+) from 134 ~ 140 5n. The E(2+) for 
140 Sn is 1949 keV showing a sudden increase for N=90, 
indicating a closed shell structure for 140 Sn. The trend 
is very similar to that observed for neutron - rich iso- 
topes of Ca while filling up the z^l/7/2) orbital and that 
shown by neutron - rich oxygen isotopes while filling up 
the ^(ld 5 / 2 ) single particle orbital (Fig.l). 

In order to put forward further evidence and to un- 
derstand the shell closure at 140 Sn more precisely, the 
effective single-particle energies (ESPE) @ for the neu- 
tron orbitals for the two Hamiltonians have been com- 
pared. The ESPE is defined as bare single particle energy 
(spe) added with the monopole part of the diagonal two 
body matrix elements (tbme) . The bare spe is originated 
from the interaction of a valence nucleon with the doubly 
closed core. The monopole interaction contribution is the 
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FIG. 2: The ESPEs of the neutron single particle orbitals for 
CWG and SMPN with increasing valence neutron number for 
Z=50. 

(2 J + 1) weighted average of the diagonal tbme, which 
arises from the interaction of a valence nucleon with the 
other valence nucleons. The ESPE for the configurations 
^(2/7/2)" in 132-140 Sn with valence neutron number n 
varying from to 8 are shown in Fig. 2. For both SMPN 
and CWG, the energy gap between v{2f 7 / 2 ) and f(3p3/ 2 ) 
single particle orbitals is 854 keV for 132 Sn core. But the 
gap between the corresponding ESPEs increases to 2.246 
MeV at N=90 with SMPN. This gap is sufficient to make 
140 Sn a doubly-magic nucleus. For CWG this gap does 
not show any increase but instead decreases slightly to 
826 keV. However shell closure at N=90 with SMPN does 
not contradict the experimentally observed fact that it is 
suitable for onset of deformation for nuclei above 132 Sn 
( with Z> 54, like Xe, Ba etc.) @ GJ, III • Fig. 3 shows 
that the proton ESPEs for SMPN favours the onset of 
collectivity at N=90 for Z>54. This is evidenced by the 
substantial reduction of the ^(157/2) and ir(2d 5 / 2 ) energy 
gap with u(2f 7 / 2 ) s . This is very similar to the appear- 
ance of the new shell gaps for the oxygen isotopes which 
disappears at larger Z values [6j. 

To analyse the origin of this new shell closure, the 
important physical aspects of both the residual interac- 
tions are extracted by a spin-tensor decomposition [2ll | of 
the two body matrix elements (tbmes). The nomencla- 
ture for the separated interaction components has been 
adopted from Ref. [22]. They are central, antisymmet- 
ric spin-orbit (ALS), spin-orbit (LS) and tensor. Fig. 4 
shows this decomposition. For SMPN, the central and 
ALS part for 2/ 7 / 2 - 2/7/2 tbmes account for majority of 
the downward shift of the ESPE of 2/7/2 with increasing 
valence neutron number (n). The tbmes involving 3^3/2 
are not modified in SMPN. They have dominant contri- 
bution from the central part. The central parts of 2f 7 / 2 
and 3p3/2 vary with similar slopes for increase in n. So 
the variation in ALS part is primarily responsible for this 
observed shell gap at N=90. 
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FIG. 3: The ESPEs of the proton single particle orbitals for 
SMPN with increasing proton number for N=90. 



Even though the neutron-neutron tbmes involving 
3p3/2 orbital can not be adjusted at present due to lack 
of data, it can be safely assumed that the ESPE of this 
orbital will not have a steeper down-sloping trend than 
that of 2/7/2 closing the presently observed 2/ 7 / 2 - 3p 3 / 2 
gap at N=90. It has been found for Ca isotopes that 
the monopole contribution of / 7 / 2 - P3/2 is positive in 
contrast to the negative values of / 7 / 2 - / 7 /2 terms [TJ . 

The ALS component in the tbmes corresponds to those 
LS-coupled matrix elements which have S^S', i.e., terms 
non-diagonal in S (spin). Thus these terms do not con- 
serve total spin of the matrix elements [23|, [24| . But the 
interactions which are parity conserving and isospin con- 
serving must also conserve the total spin. Bare nucleon- 
nucleon force contains no ALS term [23|, [24| . But ef- 
fective interaction is not simply related to bare nucleon- 
nucleon force. Core polarisation corrections to the G- 
matrix give rise to non-zero but small ALS matrix ele- 
ments. A characteristic feature common to many empiri- 
cal effective interactions is strong ALS components in the 
tbmes f23|. It usually arises from inadequate constraint 
by the data. It indicates the important contributions 
from higher order renormalisation or many body forces 
to the effective interactions. In empirical SMPN such 
many - body effects might have been included in some 
way through the modification of important tbmes. 

At this point it is important to note that SM calcula- 
tions using two-body realistic interactions derived from 
the free nucleon-nucleon force fail to reproduce some shell 
closures [l| . It is now rather well established that increase 
of the ld 5 / 2 - 2si/ 2 gap for Z=8 and l/ 7/2 - 2p 3/2 gap 
for Z=20 (as a function of neutron number) , required to 
explain empirical data are not obtained in the calcula- 
tions with these interactions. It has been shown that 
the three-body forces have to be taken into account to 
reproduce these shell gaps [l], l25l - [27j . Thus many of the 
previously observed discrepancies are now solved. Otsuka 
et al. [27} have proposed a three-body delta-hole mech- 
anism to explain these shell gaps and they have shown 
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FIG. 4: The spin - tensor decomposition of neutron ESPEs 

with increasing valence neutron number. 



that three-body forces are necessary to explain why the 
doubly-magic 24 O nucleus is the heaviest oxygen isotope. 
Zuker 26] showed earlier that a very simple three-body 
monopole term can solve practically all the spectroscopic 
problems in the p, sd, and pf shells those were hitherto 
assumed to need drastic revisions of the realistic two- 
body potentials. 

As a next step therefore, in this attempt of 
analysing the new shell gap, we have incorporated a 
simple three-body monopole term in CWG as pre- 
scribed in Ref . pa. 26]. We have incorporated cor- 
rections in 2/7/2-2/7/2 and 2/ 7 / 2 -3j>3/2 tbmes simi- 
lar to those in KB3 for I/7/2-I/7/2 an d 1/7/2-2^3/2 
tbmes. The correction terms included in the tbmes 
are Vf^J 1 (CWG3)= V^Jj 1 (CWG)-110 kcV, for J=0,4 

and 6; V f J f/f =1 (CWG3)= V f J = 2 f ' T=1 (CWG) - 310 keV 

and F / ^J= 1 (CWG3)= y/;J, =1 (CWG)+300 kcV for 3=2, 
3, 4 and 5. Here f stands for 2/ 7 / 2 and r stands for 3p 3 / 2 . 
The effective single particle energies after this correction 
are plotted in Fig. 5. A stronger 2/7/2 — 3p 3 / 2 shell gap 
now appears with CWG3 compared to that with SMPN. 

With this new interaction, CWG3, the binding energy 
of 134 Sn and E(2+) energies of 134 - 140 Sn are -6.307, 
0.627, 0.562, 0.521 and 2.532 MeV, respectively. This 
trend is very similar to that obtained with SMPN, with 
an indication of a stronger sh ell g ap at N=90. Com- 
parison of the wave functions [28| of the 0+ s and 2+ 
states with CWG and SMPN showed that whereas CWG 
favours large configuration mixing conserving seniority 
as far as possible, the SMPN on the other hand favours 
purer structure of the low-lying states, characteristics of 
K2/7/2) 8 and K2/7/ 2 ) 7 (3p 3 /2) multiplets. For CWG3, 
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FIG. 5: The neutron ESPEs for CWG3 with increasing va- 
lence neutron number, (see text for details). 

the wave function composition for the 0+ s is (81.3%) 
from the ^(2/ 7 / 2 ) 8 partition, similar to SMPN. But due 
to overestimation of the up-sloping trend of f(3p 3 / 2 ) 
ESPE (Fig.4), for 2+ state, 38.7% originates from the 
v(2f 7/2 f(lh 9/2 ) 2 and 12% from K2/ 7/2 ) 6 (2/ 5/2 ) 2 . 

In conclusion we find that comparison with the sys- 
tematics of 2^~ energies of other n-rich domains and the 
spin-tensor decomposition of the two interactions estab- 



lish the new shell closure at 140 Sn. The ALS term in- 
corporates in it the contributions of many body forces 
in the empirical interaction. A large contribution of 
this term in the ESPE of 2f 7 / 2 in empirical interaction 
SMPN has been observed. It is found to be responsi- 
ble for the gap observed in SMPN results. The CWG 
indicates a weak shell closure at 140 Sn. A simple three- 
body monopole term has been included in CWG to get 
CWG3. The new CWG3 predictions showed good agree- 
ment with that from SMPN, indicating a stronger shell 
gap at N=90 for Sn isotopes as well as decreasing 2^ en- 
ergies for 134-138 Sn. This also shows, similar to that in 
sd and fp shells, three body effect plays an important 
role for shell evolution in neutron rich Sn isotopes above 
132 Sn. The anomalously depressed 2+ states in Sn iso- 
topes having N=84-88, and the new magic number for 
N=90, might have interesting consequences for the r - 
process nucleosynthesis. 
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